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La presente invention conceme une methode pour contraindre des modeles 
stochastiques representant des zones souterraines heterogenes telles que des reservoirs 
petroliers a des donnees dites dynamiques car elles varient avec les deplacements de fluide. 
Ces donnees sont par exemple des donnees de production ou encore des donnees de 
15 pression obtenues a partir d'essais de puits. 

Etat de la technique 

L'utilisation de modules stochastiques de type gaussien pour representer 
I'heterogeneite des structures souterraines est decrite par exemple par : 

- Joumel, A.G., and Huijbregts, ChJ. ; " Mining geostatistics ", Academic Press, 1978 ; 
20 ou 

- Chiles, J.P. and. Delfiner, P. : "Geostatictics - Modeling Spatial Uncertainty",' Wiley- 
Interscience Publishers, John Wiley and Sons, 1999. 

Un modele numerique de reservoir peut etre fonme a paitir d'un ensemble de 
mailles auxquelles on attribue les valeurs d'une realisation d'un modele stochastique de 
25 type gaussien ou apparente au type gaussien. Ces valeurs peuvent etre assimilees a des 
porosites ou des perm^abilites. 

La mise en conformite du modele numerique de reservoir avec les donnees 
dynamiques mesurees sur le terrain peut etre abordee sous la forme d'un probleme 
d' optimisation. Une fonction objectif ou cout, definie au prealable, quantifie I'ecart entre 
30 les donnees dynamiques mesurees pour le milieu reel et les reponses correspondantes du 
modele numerique de reservoir. Ces reponses sont calculees a I'aide d'un simulateur 
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d'ecoulement numerique. L'objectif du probleme d' optimisation est de modifier le modele 
de reservoir ou plutot la realisation associee pour minimiser la fonction objectif. Ce 
processus est iteratif ; chaque iteration implique une simulation directe des ecoulements. 
Une bonne methode d'optimisation doit permettre : a) de modifier des realisations 
5 discretisees sur un tres grand nombre de mailles ; b) d'effectuer les modifications tout en 
respectant le modele stochastique, c'est a dire que la realisation modifiee doit etre 
coherente avec le modele stochastique ; et c) de limiter le nombre de simulations directes 
d'ecoulement, car elles requierent un temps calcul considerable. 

Parmi les techniques d'optimisation connues, on peut citer le recuit simule. Cette 
10 approche est decrite, entre autres, par : 

- Gupta, A.D. et al. : "Detailed characterization of fractured limestone formation using 
stochastic inverse approaches", SPE Ninth Symposium, 1994. 

Cette technique repose sur des echanges de valeurs de la realisation entre mailles. A 
chaque echange, la fonction objectif doit etre calculee et done une simulation directe 
15 d'ecoulement executee. Ce processus requiert un nombre excessif d'iterations. En outre, 
afin de preserver la concordance entre la realisation et le modele stochastique, un terme 
supplementaire concemant le variogramme est introduit dans la fonction objectif, ce qui 
rend 1' optimisation plus delicate. 

D'autres techniques d'optimisation, plus frequemment appliquees, s'appuient sur le 
20 calcul de gradients. Plusieurs approches s'appuyant sur les gradients sont presentees par : 

- Tarantola, A., "Inverse problem theory - Methods for data fitting and model parameter 
estimation": Elsevier Science Publishers, 1987. 

Elles necessitent le calcul des gradients de la fonction objectif par rapport aux 
parametres du probleme qui sont les valeurs de la realisation en chaque maille. Les 

25 realisations sont alors modifiees en fonction de ces gradients de sorte que la fonction 
objectif diminue. Le probleme lie au conditionnement d'un modele de reservoir a des 
donnees de production n'est pas lineaire : les techniques de minimisation par des calculs de 
gradients s'emploient it6rativeraent. Apres chaque modification de la realisation, une 
simulation directe d'ecoulement est effectuee et les gradients sont recalcules. Appliquees 

3D brutalement, les methodes des gradients conduisent a un calage des donnees dynamiques. 
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mais detruisent la coherence entre le modele stochastique et la realisation. En outre, elles 
ne permettent pas de considerer un tres grand nombre de parametres. Pour eviter ces 
limites, des techniques de parametrage geostatistique peuvent y etre integrees. On peut 
citer a ce stade la m^thode des points pilotes, decrite par : 

5 - de Marsily, G. et al. : "Interpretation of interference tests in a well field using 
geostatistical techniques to fit the permeability distribution in a reservoir model" in 
Verly, G. et al. (ed.), Geostatistics for natural resources characterization. Part 2, D. 
ReidelPub. Co, 1984. 

Cette methode consiste a selectionner dans la realisation un certain nombre de 
10 points, dits pilotes, a calculer les derivees de la fonction objectif par rapport aux valeurs en 
ces points, ^ modifier les valeurs de ces points en consequence et a propager la 
perturbation ainsi definie par une technique de krigeage. La methode des points pilotes 
peut induire des variations aberrantes des valeurs des points pilotes. 

Une autre technique de parametrage geostatistique, qui evite I'ecueil precedemment 
15 mentionne, est la methode de deformation graduelle d'un modele stochastique d'un milieu 
heterogene tel qu'une zone souterraine. Elle est exposee et exploitee par Hu, L.-Y. et 
al.,dans les brevets FR 2.780.798 et FR2.795.841 du demandeur. 

La methode de deformation graduelle permet de modifier graduellement une 
realisation d'un modele stochastique de type gaussien ou de type apparente au modele 
gaussien tout en respectant ce modele. La realisation deforraee est encore une realisation 
du modele stochastique. Lorsque la methode de deformation graduelle est introduite dans 
un processus d'optimisation, la demarche est la suivante. La realisation initiale est 
combinee avec un nombre fixe de reaHsations independantes se rattachant au meme 
modele stochastique. Ces realisations sont appelees realisations complementaires. La 
combinaison est controlee par autant de parametres de deformation qu'il y a de realisations 
complementaires. Elle foumit une nouvelle realisation. Puis sont calculees les derivees de 
la fonction objectif par rapport aux parametres de deformation. Ces demiers sont modifies 
afin de prendre en compte I'information issue des derivees. Une premiere optimisation par 
rapport aux parametres de deformation foumit une realisation verifiant le modele 
stochastique el reduisant la fonction objectif En general, ce processus d'optimisation doit 
etre repete plusieurs fois avec des realisations complementaires differentes afin de 
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dimitiuer suffisamment la fonction objectif, ce qui peut signifier dans certains cas 
I'execution d'un nombre prohibitif de simulations directes d'ecoulement. 

La methode de deformation graduelle permet, dans certains cas, de modifier une 
realisation localement. Cette possibilite est fondee lorsque le deformation graduelle est 
5 combinee au generateur geostatistique FFTMA decrit par : 

- Le Ravalec, M. et al. : The FFT moving average (FFT-MA) generator: an efficient 
numerical method for generating and conditioning Gaussian simulations. Math. Geol., 
32(6), 2000. 

Ce generateur produit des realisations pour un modele stochaslique de type 
10 gaussien specific au prealable en convoluant un bruit blanc gaussien avec un operateur 
dependant de la fonction de covariance. Une deformation locale peut etre effectuee en 
appliquant la methode de deformation graduelle au bruit blanc gaussien sous-jacent a la 
realisation. 

Les techniques de gradients developpees & ce jour pour le calage par rapport aux 
15 donnees dynamiques reposent sur un lien direct entre les variations auxquelles est soumise 
la realisation representant le reservoir et la variation de la fonction objectif. Le processus 
d' optimisation implique de modifier d'abord la realisation et de lancer ensuite une 
simulation d'ecoulement pour apprehender les variations qui en resultent pour le calage. 

Une approche differente, valable dans le cas ou les ecoulements sont modelises par 
20 des lignes de courant, est proposee par : 

Wang, Y., and Kovscek, A.R. : A streamline approach for history-matching production 
data. SPE/DOE lOR, 2000. 

Inversant la demarche des approches traditionnelles, ces auteurs proposent une 
m^thodologie fondee sur Tapplication successive de deux etapes. La premiere etape se 

25 concentre sur le calage : elle permet d'evaluer la modification a apporter aux permeabilites 
effectives des lignes de courant pour ameliorer le calage. La deuxieme etape conceme le 
transfert de la variation de permeabilite effective des Ugnes de courant vers la realisation. 
Ce processus est poursuivi jusqu'a Tobtention d'un calage juge satisfaisant. La limite 
majeure de cette approche est de ne pas preserver la coherence de la realisation vis a vis du 

30 modele stochastique. 
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La methode selon I'invention 

La methode selon I'invention permet d'elaborer des modeles numeriques, 
representant des milieux heterogenes souterrains tels que des reservoirs petroliers ou des 
aquiferes, conformes tout d'abord avec un ensemble de donnees dynamiques mesurees 
5 dans des puits de production, d'injection ou d'observation, mais aussi avec un modele 
stochastique qui caracterise la distribution spatiale des h^terogeneites. 

Elle comporte essentiellement una minimisation de la fonction objectif a partir d'un 
processus iteratif subdivise en deux etapes et par I'integration de la methode de 
deformation graduelle dans la deuxieme etape. La premiere etape identifie la perturbation h. 

10 apporter aux permeabilites effectives de zones predefinies pour reduire la fonction objectif. 
La deuxieme etape se concentre sur la minimisation de I'ecart entre les permeabilites 
effectives souhaitees et les permeabilites effectives correspondantes calculees pour la 
realisation consideree, en modifiant cette meme realisation a I'aide de la methode de 
deformation graduelle. Cette etape de minimisation est une optimisation intermediaire, qui 

15 ne requiert aucune nouvelle simulation d'ecoulement et qui permet de deformer la 
realisation en respectant le modele stochastique gaussien. Si la fonction objectif n'est pas 
satisfaisante a ce stade, le processus reprend au niveau de la premiere etape. 

La methode selon I'invention est une methode d'optimisation graduelle decouplee. 
Le terme graduel s'applique dans la mesure ou la deformation du modele numerique se fait 

20 suivant les principes fondamentaux de la methode de deformation graduelle. Le terme 
decouple rappelle que 1' optimisation s'appuie sur deux etapes, la premiere definissant la 
correction a apporter aux permeabilites effectives, la deuxieme la propageant 
independamment vers le modele numerique. Ces deux etapes impliquent que les 
simulations d'ecoulement sont executees uriiquement pendant la premiere etape afin 

25 d'ameliorer le calage des permeabilites effectives. En aucun cas, la deuxieme etape ne 
requiert de simulation d'ecoulement pour estimer la maniere dont doit etre deforme le 
modele numerique de reservoir. 

En d'autres termes, la methode selon I'invention permet de former iterativement un 
modele representadf du champ de permeabilite d'un milieu heterog^ne, discretise par un 
30 maillage, tel qu'une zone souterraine, contraint par des donnees geologiques a priori et des 
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donnees dynamiques collectees dans le milieu par des mesures et observations prealables. 
Elle comporte: 

- une premiere etape dans laquelle: 

a) on genere un modele initial de la permeabilite en adequation avec un modele 
5 stochastique gaussien ou apparente, coherent avec les donnees geologiques a priori, et on 

realise au moyen d'un simulateur appropri6, une simulation des ecoulements de fluides ; et 

b) on identifie des zones a I'interieur du reservoir, on calcule les permeabilites 

effectives de ces zones, et, a partir des resultats du simulateur, on estime les corrections a 
apporter a ces permeabilites effectives pour ameliorer le calage par rapport aux donnees; et 

- une deuxieme etape dans laquelle c) on propage les corrections a I'ensemble des 
mailles du champ de permeabilite, par un processus d'optimisation iteratif dans lequel on 
minimise une fonctionnelle qui quantifie la difference entre les permeabilites effectives 
necessaires pour obtenir le dit calage et les permeabilites effectives calculees pour le 
champ de permeabilite considere, en utilisant une technique de deformation graduelle de 
realisations du modele stochastique. 

On definit les zones soit manuellement, soit automatiquement a partir du simulateur 
d'ecoulement. 

Suivant un mode de mise en ceuvre, on realise la simulation des ecoulements au 
moyen d'un simulateur par lignes de courant, et on identifie les zones du milieu par un 
20 ensemble de mailles traversees par une ou plusieurs lignes de courant, de geometrie fixe. 

Suivant un autre mode de mise en ceuvre, on identifie les dites zones sous la forme 
de portions de volume a la peripheric de puits traversant le dit milieu, dans le cadre 
d'essais de puits. 

Suivant un autre mode de mise en ceuvre, on attribue au moins un parametre de 
25 deformation graduelle a chacune des zones. 

Par rapport a I'art anterieur, la methode selon I'invention permet de caler une 
realisation d'un modele stochastique a des donnees dynamiques en gardant la coherence 
par rapport au modele stochastique et en reduisant significativement le nombre de 
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simulations directes d'^coulement a executer. Contrairement aux approches developpees a 
ce jour, urie deformation de la realisation n'implique pas systematiquement une nouvelle 
simulation d'ecoulement. Elle donne ainsi la possibilite de mieux explorer I'espace des 
realisations et de determiner rapidement non pas une realisation contrainte, mais plusieurs. 
5 De plus, la methode selon I'invention permet de deformer la realisation a partir d'un grand 
nombre de parametres ce qui acc61^re Toptimisation intermediaire. Elle rend aussi possible 
la deformation de la realisation par zones, celles-ci pouvant 8tre les memes que les zones 
predefinies pour le calcul des permeabilites effectives lors de la premiere etape. 

Presentation des figures 

10 D'autres caracteristiques et avantages de la methode selon I'invention apparaTtront 

plus precisement a la lecture de la description ci-apres d'un exemple non liinitatif 
d'application, en se referant aux dessins annexes oia : 

- la Fig. 1 montre un exemple de schema d'optimisation en deux etapes selon la 
methode ; et 

15 - la Fig. 2 montre une comparaison entre un flux fractionnaire reel au flux fractionnaire 
simule correspondant. 

Description detaillee de la methode 

Traditionnellement, le calage des modeles de reservoir est un processus iteratif au 
cours duquel on execute une simulation d'ecoulement a chaque fois qu'on perturbe le 

20 reservoir. Par souci de simplification, on considere que le modele numerique de reservoir 
se ramene a une realisation d'un modele stochastique pour la permeabilite. En d'autres 
termes, une simulation d'ecoulement est requise pour toule variation du champ de 
permeabilite. Au contraire, la methode selon I'invention peut proposer plusieurs variations 
du champ de permeabilite en ne recourant qu'^ une unique simulation d'ecoulement. Pour 

25 ce faire, elle procede en deux etapes (Fig. 1). Lors de la premiere etape, on determine des 
zones du rdservoir et on estime la correction a apporter aux permeabilites effectives de ces 
zones pour ameliorer le calage des donnees. Puis, un processus d'optimisation est mis en 
place pour propager la perturbation determinee pour les permeabilites effectives de ces 
zones vers le champ de permeabilite representant le reservoir. La premiere etape requiert 

30 une simulation d'ecoulement puisqu'elle depend d'une comparaison des donnees reelles 
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avec les reponses synthetiques correspondantes. Au contraire, la deuxieme etape ne 
demande pas de simulation d'ecoulement supplementaire. Enfin, 1' optimisation qui 
intervient dans la deuxieme etape est parametree suivant la methode de deformation 
graduelle ce qui permet de preserver la coherence de la distribution des valeurs de 
5 permeabilite (ou toule autre propriety consideree) vis a vis du modele de variabilite 
spatiale. 

L'algorithme general est resume comme suit : 

Un champ initial de permeabilite en adequation avec un modele stochastique ayant 
ete genere ; 

10 a) etape 1 - on fait un simulation d'ecoulement, on calcule les permeabilites 

effectives de ces zones, on identifie des zones a I'interieur du champ considere et on estime 
la variation a appliquer aux permeabilites effectives de ces zones afin de reduire I'ecart 
entre les donnees de production mesurees sur le terrain et les donnees de production 

simulees ; 

15 b) etape 2 - la perturbation requise au niveau des permeabilites effectives des 

zones est propagee a I'ensemble du champ de pemneabilite par un processus 
d'optimisation. La fonction a minimiser quantifie la difference entre les permeabilites 
effectives souhaitees et les permeabilites effectives calculees pour le champ de 
permeabilite considere. La technique de deformation graduelle est mise en oeuvre pour 

20 modifier le champ de permeabilite. 

c) retourenb)tantquelecalagen'estpassatisfaisant. 

Le calcul des permeabilites effectives des zones peut etre fait en utilisant le 
simulateur ou bien encore etre obtenu par utilisation de techniques de mioyenne bien 
connues des gens de I'art. 

25 La methode selon I'invention est adaptable a tout simulateur d'ecoulement dans la 

mesure ou Ton se donne les moyens de definir des permeabilites effectives pour differentes 
zones du reservoir, ces demieres etant identifiees manuellement par I'utilisateur ou 
automatiquement a partir d'un critere donne. Par exemple, pour un simulateur de tests de 
puits, on pourra selectionner des couronnes de rayon croissant, centrees sur les puits, pour 
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definir les zones. Les penneabilites effectives pourront dans ce cas etre reliees aux 
permeabilites apparentes. 

Dans I'exemple developpe a present, nous nous concentrons sur le cas d'un 
simulateur par lignes de courant. Les zones considerees dans ce cadre sont ddterminees par 
5 les lignes de courant elles-memes. Chacune des deux etapes precedemment mentionnees 
est revue dans le detail dans les deux sections qui suivent. 



L'objectif de la premiere partie est de decrire la geometric des lignes de courant et 
d'identifier la modification a apporter a la permeabilite effective de ces lignes pour mieux 
10 caler les donnees de production. Une ligne ou un ensemble de lignes de courant definissent 



La modelisation des ecoulements par des lignes de courant presentent certaines 
caracteristiques. Le deplacement de fluide le long d'une ligne de courant est un probleme 
unidimensionnel ; les lignes de courant ne communiquent pas entre elles. Quand le rapport 
15 des mobilites des fluides en presence est unitaire et quand les conditions aux limites ne 

changent pas, la georaetrie des lignes de courant est fixe. Quand le rapport des mobilites 
est different de un, on distingue deux alternatives : on fixe la geometrie des lignes de 
courant et on autorise les variations des flux pendant le deplacement des fluides ou on 
distribue le flux uniformement entre les tubes de courant et on remet a jour periodiquenient 
20 la geometrie des lignes de courant. Nous nous plagons d'emblee dans la premiere 
configuration. 

Pour simplifier le probleme, on considere des groupes de lignes de courant plutot 
que les lignes de courant individuellement : on reduit ainsi le nombre de parametres. La 
permeabilite effective d'un groupe de lignes de courant s'exprime comme la moyenne 
25 harmonique ponderee par les volumes de fluide des permeabilites des mailles traversees 



Etape 1 : Modification des permeabilites effectives 



une zone. 



par les lignes de courant : Kf^ = 



Eq. 1 
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xf-^ est la permeabilite effective du i^""^ groupe de lignes de courant. N^sl est le 
nombre de lignes de courant comprises dans ce i^""^ groupe, Nk est le nombre de mailles 
traversees par la k^""" ligne de courant du i^""^ groupe. q,^ j est le flux pour la k^™" ligne de 
courant au niveau de la j^™*" maille. Ar^ y est le temps de vol pour la ligne de courant 
5 au travers de la j"^""* maille. 

Pour une realisation de depart donnee, la simulation d'ecoulement pemnet de 
determiner la geometric des lignes de courant et de calculer leurs peimeabilites effectives. 
II reste a evaluer la variation de permeabilite effective qui foumirait une amelioration du 
calage des donnees de production. Pour ce faire, on compare, aux puits producteurs, les 
flux fractionnaires observes et les flux fractionnaires correspondants simules (Fig. 2). On 
ordonne alors les lignes de courant dans I'ordre croissant de leurs temps de percee et on 
discretise les courbes de flux fractionnaires. Les segments ainsi definis sont associes a 
des groupes de lignes de courant. En I'absence de concordance entre les flux simules et les 
flux reels, les permeabilites effectives des groupes de lignes de courant sont considerees 
comme respon sables des differences. Une coirection appliquee a 1' ensemble des groupes 
de lignes de courant permet de restreindre cet ecart. Considerons le groupe de lignes de 
courant associes a un increment de flux A(y,- (Fig. 2). La permeabilite effective qu'on 
souhaiterait avoir pour le groupe i de lignes de courant afin d'ameliorer le calage est : 

j^eff _ h.simule ^eff 

i,souhaite i.simule 

20 ou ti,reference ©t ti_simuie sont Ics temps de pcrcec moycns associes a I'increment A^; pour les 
courbes de flux fractionnaires reel et simule. On precede de maniere analogue pour chaque 
groupe de lignes de courant. 

Etape 2 : Optimisation intermediaire - Modification du champ de 
permeabiiite 

25 La premiere etape determine une variation a appliquer aux permeabilites effectives 

des lignes de courant pour caler les flux fractionnaires. L'objectif de la deuxieme etape est 
de transferer cette variation des lignes de courant vers ie champ de permeabilite tout en 
preservant la coherence du champ de permeabilite vis a vis du modele stochastique. Pour 
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ce faire, on definit un probleme d'optimisation. II s'agit ici de minimiser la fonction, dite 
intermediaire : 

Ng est le nombre de groupes de lignes de courant. Kff^^^^-^^&st la permeabilite 

5 effective souhaitee pour le groupe i de lignes de courant afin de reduire les differences 
entre les flux fractionnaires mesures et simules : cette valeur a ete determinee lors de la 
premiere etape. Kf^-^^^i^ est la permeabilite effective du groupe / de lignes de courant pour 

le champ de permeabilite considere. Ce probleme d'optimisation est non lineaire et peut 
impliquer autant de parametres que le champ de permeabilite comprend de mailles. En y 

10 integrant la methode de deformation comme technique de parametrage, on se donne la 
possibilite de reduire le nombre de parametres et on assure une modification conforme du 
champ de permeabilite par rapport au modele stochastique definissant la distribution 
spatiale des heterog^neites dans le champ de permeabilite. De plus, afin de pouvoir 
modifier le champ de permeabilite par zones, les zones correspondant aux differents 

15 groupes de lignes de courant, la deformation graduelle est appliquee non pas au champ de 
permeabilite lui-meme, mais au bruit blanc gaussien sous-jacent (Fig. 1). 

La version la plus simple de la methode par defonnation graduelle consiste a 
combiner deux bruits blancs gaussiens, Z; et Z2, suivant la relation : 

z(x)(p) = Zi (x)cos(?rp) + ^2(x)sin(7r/?) , 
20 oil p est le parametre de deformation et x la position. Le gradient de z par rapport a /? est : 

^ = -TTZi sin(;ryo) + 71Z2 cosinp) . 
dp 

Le bruit blanc gaussien z est ensuite transforme en une realisation gaussienne y de 
modele de covariance C, de moyenne nulle et variance unitaire, a I'aide d'un produit de 
convolution, developpe au coiur du generateur FFTMA : 
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/ resulte de la decomposition de la fonction de covariance. La derivee de y par 
rapport au parametre de deformation est : 

dp dp' 

Tout autre generateur geostatistique produisant des realisations gaussiennes ou de 
5 type apparente peut 6tre employe a la place du generateur FFTMA ^ condition d'y integrer 
le calcul des gradients. Le generateur FFTMA off re neanraoins I'avantage de la rapidite 
d'execution, meme pour des realisations discretisees sur un tres grand nombre de mailles. 

La realisation y, pour I'instant centree reduite, c'est a dire de moyenne nuUe et 
variance unitaire, peut ensuite etre transformee en une realisation gaussienne w de 
10 moyenne m et variance (T : 

w(x) = m + ayix) 
La derivee par rapport au parametre de deformation devient : 

A ce stade, il faut aussi prendre en compte des donnees staliques qui auront ete 
15 observees en des points precis, par exemple dans des puits. Cette infomiation est en 
general integree dans la realisation generee a I'aide d'une technique de krigeage. La 
realisation contrainte Wc se deduit de : 

(x) = (x) + (vw(x) - (xj) 

ou WdK et Vi'K sont les realisations obtenues pour la premiere a partir du krigeage des 
20 observations reelles et pour la seconde du krigeage des valeurs de w aux points 
d' observation. L'estimateur de krigeage, dans le cadre dual, s'exprime : 

WifW = Sp.C'(x-x,) + m 



C est la fonction de covariance. Les x, sont les positions des n observations. Les 
poids Pi sont independants de la position, mais dependent des parametres de deformation. 
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On peut montrer que la derivee de w par rapport au parametre de defomiation s'obtient a 
partir : 



Le gradient des permeabilite par rapport au parametre de deformation s'exprime : 



L'equation 1) permet de calculer la permeabilite effective d'un groupe de lignes de 
courant pour un champ de permeabilite donne. On en deduit la derivee de la permeabilite 
effective par rapport au parametre de deformation : 



Ces differentes relations montrent comment deformer un champ de permeabilite et 
calculer les gradients des permeabilites effectives des lignes de courant par rapport a un 
parametre de defoirnation. Toutes ces relations sont aisement generalisables au cas ou plus 
parametres de deformation interviennent. On poun-a notamment faire le choix d'affecter un 
parametre de deformation a chacune des zones definies pour le calcul des permeabilites 
effectives. 




Le champ de permeabilite k se deduit de la transformation log-normale de Wc : 



A:(x) = exp(w^(x)). 



dp^ ^~ dp 



■k{x) 




Plusieurs techniques peuvent eire envisagees pour minimiser la fonction 
intermediaire, mais puisque les gradients des permeabilites effectives sont disponibles pour 
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un cout calcul negligeable, une approche de type Gauss-Newton peut s'averer judicieuse. 
Pour determiner la perturbation a apporter aux parametres de deformation pour reduire la 
fonction intermediaire, resout le systeme : 

HAp = -V(F/) 

Ap est la perturbation a definir, V(F/) comprend les gradients de la fonction 
intennediaire par rapport aux parametres de defonnation et H est une matrice approchee de 
la matrice hessienne : 

H = G'WG 

G est la matrice de sensibilite : elle rassemble les derivees des permeabilites 
effectives des groupes de lignes de courant par rapport aux parametres de deformation. W 
est la matrice poids : elle egale ici la matrice identite. 

Finalement, I'algorithme developpe pour minimiser la fonction intermediaire est 

decrit comme suit. L'etape 1 permet de definir d'une part les permeabilites effectives 
souhaitees pour les groupes de lignes de courant et d' autre part des zones affectees a ces 
groupes. 

a) On affecte au moins un parametre de deformation par zone, 

b) On definit un bruit blanc gaussien initial et au moins un bruit blanc gaussien 
complementaire. 

c) On effectue une deformation graduelle en combinant les bruits blancs 
gaussiens precedents. On calcule dans le meme temps les derivees du bruit blanc gaussien 
z issu de la defonnation graduelle par rapport aux parametres de deformation. 

d) On transforme le bruit blanc gaussien z en un champ de permeabilite k, on 
detennine aussi les derivees des permeabilites par rapport aux parametres de deformation. 

e) On calcule les permeabilites effectives des groupes de lignes de courant 
definies lors de l'etape 1, ainsi que leurs derivees. 

f) On estime alors la perturbation a apporter aux parametres de deformation 
pour reduire la fonction intermediaire. 



g) A ce stade, plusieurs alternatives se dessinent. Si la fonction intermediaire 
n'est pas assez faible et si elle ne semble pas avoir converge, on remet a jour les 
parametres de deformation et on revient en c). Si la fonction intermediaire n'est pas assez 
faible, mais semble avoir converge, on revient en b), c'est a dire qu'on met a jour le bruit 
5 blanc gaussien initial et qu'on tire au hasard un nouveau bruit blanc gaussien 
complementaire. Enfin, si la fonction objectif est suffisamment faible ou si elle ne semble 
pas vouloir decroitre davantage, I'etape 2 prend fin. 

La procedure decrite ici se concentre sur le cas des simulations d'ecoulement par 
ligne de courant. hes lignes de courant apparaissent en effet comme un outil tres naturel 

10 pour definir des zones. Le coeur de la presente invention impliquant I'identification de 
zones et le calcul des permeabilites effectives pour ces zones, ce choix semble logique. 
D'autres types d' applications sont neanmoins tout a fait envisageables. Par exemple, 
rinfomiation relative aux differentes zones poun'aient etre foumie par des tests de puits. 
La permeabilite effective peut dans ce cas etre calculee en fonction du rayon 

15 d' investigation autour du puits ; les zones considerees sont alors des couronnes. Le 
simulateur d'ecoulement mis en oeuvre pour simuler les essais de puits peut etre un 
simulateur d'ecoulement standard. La procedure a suivre est analogue a celle decrite pour 
les lignes de courant. La simulation d'ecoulement permet d' identifier des zones et de 
determiner la permeabilite effective pour ces zones qui peut etre comparee aux donnees 

20 mesurees sur le terrain. Puis, la minimisation d'une fonction objectif intermediaire suivant 
I'approche presentee plus haut permet de propager la correction a apporter aux 
permeabilites effectives vers les permeabilites absolues des mailles a I'interieur des zones 
tout en respectant le modele de variabilite spatiale a priori. 



REVENDICATIONS 



1) Methode pour former iterativement un modele representatif du champ de 
permeabilite d'un milieu heterogene, discretise par un maillage, tel qu'une zone 
5 souterraine, contraint par des donnees geologiques a priori et des donnees dynamiques 
collectees dans le milieu par des mesures et observations prealables, caracterisee en ce 
qu'elle comporte: 

- une premiere etape dans laquelle: 

a) on genere un modele initial de la permeabilite en adequation avec un modele 
10 stochastique gaussien ou apparente, coherent avec les donnees geologiques a priori, et on 

realise au moyen d'un simulateur approprie, une simulation des ecoulements de fluides ; et 

b) on identifie des zones a I'interieur du reservoir, on calcule les permeabilites 
effectives de ces zones, et, a partir des resultats du simulateur, on estime les corrections a 
apporter a ces permeabilites effectives pour ameliorer le calage par rapport aux dites 

15 donnees; et 

- une deuxieme etape dans laquelle c) on propage les dites corrections a I'ensemble 
des mailies du champ de permeabilite, par un processus d'optimisation iteratif dans lequel 
on minimise une fonctionnelle qui quantifie la difference entre les permeabilites effectives 
necessaires pour obtenir le dit calage et les permeabilites effectives calculees pour le 

20 champ de permeabilite considere, en utilisant une technique de deformation graduelle de 
realisations du modele stochastique. 

2) Methode selon la revendication 1, caracterisee en ce que les zones sont definies 
soit manuellement, soit automatiquement a partir du simulateur d'ecoulement. 

3) Methode selon la revendication 1, caracterisee en ce que I'on realise la 
25 simulation des ecoulements au moyen d'un simulateur par lignes de courant, et on identifie 

les zones du dit milieu par un ensemble de mailies traversees par une ou plusieurs lignes de 
courant, de geometric fixe. 
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4) Methode selon la revendication 1, caracterisee en ce que Ton identifie les dites 
zones sous la forme de portions de volume a la peripherie de puits traversant le dit milieu, 
dans le cadre d'essais de puits. 

5) Methode selon Tune des revendications precedentes, caracterisee en ce que Ton 
attribue au moins un parametre de deformation graduelle a chacune des dites zones. 
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